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Abstract 



We study the large-time behavior of solutions to the heat equation in nonhomogeneous media 
with critical singular density 

CM ' 

\x\~ 2 d t u = Au, in R N x (0,oo) 

'"^ ■ in dimensions N > 3. The asymptotic behavior proves to have some interesting and quite striking 

properties. We show that there are two completely different asymptotic profiles depending on 
whether the initial data uq vanishes at x = or not. Moreover, in the former the results are true 
only for radially symmetric solutions, and we provide counterexamples to convergence to symmetric 
■ profiles in the general case. 
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1 Introduction 

The aim of this work is to establish the asymptotic behavior of solutions to the following heat equation 
in nonhomogeneous media with critical density: 

\x\- 2 d t u{x,t) = Au(x,t), (x, t) G K N x (0, oo). (1.1) 

as a part of an ongoing project of studying the asymptotic behavior for 

\x\- 2 d t u(x, t) = Au m (x, t), (x, t)£R N x (0, oo), (1.2) 
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with m > 1. The technically more involved problem with the large-time behavior for m > 1 will be 
studied in a forthcoming paper [8]. 

Equations of type (|1.2p with general densities have been proposed by Kamin and Rosenau in a 
series of papers [llj \12\ \13\ to model thermal propagation by radiation in non-homogeneous plasma. 
Since then, many papers were devoted to developing rigorously the qualitative theory or asymptotic 
behavior for 

g(x)d t u(x, t) = Au m (x, t), (x, t) G R N x (0, oo), (1.3) 

usually asking that g{x) ~ as \x\ — > oo, for some 7 > (e.g. [18]). Thus, for m > 1, it has 

been noticed that while < 7 < 2, the solutions have similar properties to the ones of the pure porous 
medium equation (for short, PME) 

d t u = Au rn , 

see [171 119j . while for 7 > 2, the properties of the solutions depart strongly from the ones of the PME 
|10| . Thus, 7 = 2 is critical, and the asymptotic behavior for this case is left open in |10j with a 
conjecture giving the explicit profile. On the other hand, the asymptotic behavior for (jl.3|) with a 
density g{x) ~ \x\~ 2 at infinity, but g regular near the origin, is studied in another recent work |15| . 
obtaining an explicit profile solving (|1.2|) . and proving the convergence towards it in the outer region 
(that is, outside small compacts near x = 0). Unusually, the linear diffusion problem m = 1 has been 
studied later than its nonlinear version, see for example [2J[3]. 

In all cases, it has been noticed that the solutions converge asymptotically towards profiles coming 
from the pure power density equation, that is, g(x) = |x|~ 7 . Moreover, 

\x\-tdtu(x,t) = Au m (x,t) (1.4) 

has some more interesting feature: a singularity at x = 0, apart from the decay at infinity. Recently, in 
[7] the authors study formally some properties of radial solutions to f)1.4[) as a first step to understand 
its general behavior. Moreover, a study of existence and uniqueness for (|1.4p for 7 > 2 and m > 1 is 
done in [TOj Section 6]. 

Coming back to our problem, that of letting m = 1, 7 = 2 in (jl.4p . we find explicit asymptotic 
profiles which explain better the effect of the singularity at x = 0. The general case m > 1 will be 
treated in the companion paper [8]. But in order to explain these comments and to make precise the 
motivation for this work, let us state the main results of the paper. 

Main results. We deal with the Cauchy problem associated to Eq. (jl.ip with initial data 

n G Ll(R N ), n > 0, (1.5) 
where the dimension is N > 3 (except when specified) and 

Ll(M N ) := \ h : R N !->• R, h measurable, / \x\~ 2 h(x) dx < 00 i . 

I Jr n J 

In Section [2] we give the precise notions of weak solution and strong solution to (jl.ip and we prove 
that the Cauchy problem for (jl.ip is well-posed in our framework. We refer the interested reader to 
Definition 12.11 and Theorem 12.21 for the precise statements. 

We state the results about the large-time behavior, that we find quite interesting and unexpected. 
We begin with the case when uo(0) = 0. 

Theorem 1.1. Let u be a radially symmetric solution of Eq. (jl.ip with initial data satisfying f 1 1 . 5 1) 
and moreover 

M uo := / \x\~ N u (x)dx < 00, n (0)=0. (1.6) 
Jr n 
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u(x, t) 



Then we have 

lim t l/2 

t— >oo 

where ui\ is the area of the unit sphere in ]R JV and 
F(x,t) : = 



Muc 



^ r log\x\+(N-2)t 
4nt V 2^ 



■F(x,t) 



G(0 



for \x\ ^ 0, 



(1.7) 



(1.8) 



0, for x = 0, 

Remarks, (a) Let us notice also that 

max{F(x,t) : x / 0} = 0(r 1/2 ), as t ->■ oo, 

showing that the time-scale i 1 / 2 is the correct one for the asymptotic behavior in (jl.7p . More precisely, 

1 



(1.9) 



as it will be analyzed in the remarks at the end of Section [H 

(b) The mass M UQ in (|1.6p is conserved along the flow. This will be obvious from the proof of Theorem 

m 



Counterexamples in the non-radial case. We emphasize on the fact that this is true only for 
radially symmetric solutions. For general solutions, Theorem II .11 is not true. It is quite surprising that, 
due to the singularity at x = 0, the solutions that start from uq(0) = do not converge asymptotically 
to a radial profile. 

We construct a counterexample to Theorem 11.11 for solutions that are not radially symmetric. We 
pass to generalized spherical coordinates x = (r, (f>i, 4>n-2, 0) in ^ N and define the following function 

F N {x,t) =F N {r, </>!,.. .,<t>N- 2 ,6) :=0F(r,t). 

Using the formula of the Laplace operator in spherical coordinates, from the particular form of Fn 
(that involves no dependence on (0i, 4>n-2) and 8 2 Fn = 0), we find that Fn is a solution to (jl.ip . 
since 

AF N (x,t) = 9A r F(r,t) = e\x\~ 2 F t {r,t) = \x\~ 2 d t F N (x,t), 

where A r is the Laplacian operator in radial variable. The fact that Theorem 11.11 is not satisfied on 
this example is obvious since (jl.7p becomes 

lim t 1 / 2 \\(9-K)F(r,t)\\ oo =0, 

t— >oo 

for some constant K depending on the dimension N, which is false since 8 is variable, K is constant 
and t 1 /2|| F (. ; i)|| oo = i/yfa. 

When the value of the initial data at the origin is nonzero, things are completely different, as the 
following theorem states. 

Theorem 1.2. Let u be a general solution of Eq. (jl.ip with initial data uq satisfying (|1.5p and 
uq(0) = K > 0. // there exist 5, e > such that the initial data uq satisfies 

\K — uq{x)\ < \x\ 5 , as \x\ — > 0, uq(x) < \x\~ 



as \x\ 



oo, 



then we have 



where 



E(x,t) :-- 



u(t) - ^E{t) 



erfc( los|a;| 2 + Jf- 2)t 
K, for x = 0. 



= 0(r 1/2 ), ast-^oo, 

DO 

, eric(0 = ^f™e- 92 de, for \x\ ^ 0, 



(1.10) 
(1.11) 

(1.12) 
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Remarks. 1. Theorem 1 1 . 2 1 can be improved by relaxing and generalizing the condition (jl . 1Q|) : indeed, 
p. lip holds true if there exists some functions <l?i, <3?2 : (0,oo) h-» (0, 00) such that 

$i(r) 

dr < 00, / - . , dr < 00, 



'0+ r J r$ 2 (r) 

and 

uq{x) < $i(|x|), as x — > 0, uo(x) < $2(1^1); as \x\ — > 00. 
The proof is totally similar. 

2. Theorem 11.21 is indeed true in general, also without radial symmetry, in contrast with Theorem 
11.11 These two different results show the importance of the singularity of the density |x|~ 2 . Indeed, the 
evolution gives rise in a striking way to two different asymptotic profiles, only depending on a difference 
in the pointwise value of uo at the origin. 

When dealing with radially symmetric solutions, the results can be improved both with respect to 
the condition on the initial data uq and with respect to the rate of convergence. This is gathered in 
the following 

Theorem 1.3. Let u be a radially symmetric solution for Eq. (jl.ip whose initial data satisfies 

(EE}. 

(a) The asymptotic convergence (jl.lip holds true under the following condition on uq: 

h ■= / \x\~ N \K -u (x)\dx + / \x\~ N \u {x)\dx < 00. (1.13) 

Jb(o,i) Jwl n \b(q,i) 

(b) If we ask furthermore the following condition on the gradient 

h:= t K b g^l) 3 l |vn (x)|(ix<oo, (1.14) 

Jr n m 

then we have a rate of convergence: 

= 0(t~ 3/2 ), ast^oo. (1.15) 



u(t) - ^E(t) 



We will add here two more interesting consequences about the behavior near x = 0, where the 
singularity at x = plays a key role. 

Proposition 1.4. (a) There is no fundamental solution of Eq. (|1.1|) . that is, a solution with initial 
data uo = MSq. 

(b) If the initial condition uq satisfies uq(x) = for any x G -B(0, tq), for some ro > 0, then u(x, t) > 
for all (x, t) G (R N \ {0}) x (0, 00) and u(0, t) = for any t > 0. 

We plot the asymptotic profiles in Figure [1] in dimension N = 3, the general picture being similar. 
In the pictures we see the profiles evolving at various times. 

Organisation of the paper. Before passing to the asymptotic behavior, we begin by studying the 
well-posedness for (jl.ip with suitable initial data. This is the goal of Section [2J which insures us that 
the object of our study exists in suitable spaces. We then prove the main results in two steps. In a 
first step, we describe a transformation mapping radially symmetric solutions of (|1.1|) into solutions of 
the one-dimensional heat equation. This is done in Section [3l With the aid of it, we prove Theorem 
11.11 and part (a) of Theorem 11.31 Then, in order to prove part (b) of Theorem 11.31 we need one 
more transformation step. The proofs of the theorems for radially symmetric solutions and some more 
remarks about the asymptotic profiles are the subject of Section |U Finally, in Section [5] we prove 
Theorem 11.21 using as essential tool a comparison principle proved in Section [21 and we end with the 
proof of Proposition 11.41 We close the paper with a section of comments where we include a brief 
discussion of the dimension N = 1 for (jl.ip and raise some open problems. 
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Figure 1: Profiles E and F for dimension N = 3 at various times. 

2 Existence and uniqueness 

Before proving the main results about asymptotic behavior, we have to develop a theory of existence 
and uniqueness for the Cauchy problem 

\x\~ 2 u t = An, in R N x (0,oo), , . 

u(x,0) = u {x), xeR N , ^ ' ' 

for suitable initial data no £ L\(R N ). We start from the precise definition of a solution. 

Definition 2.1. We say that a function u(x,t) is a weak solution to (|2.ip if it satisfies the following 
conditions: 

(a) u e C([0,co);Ll(R N ))nL oo (R N x (r,oo)), for any t > 0, and u > inR N x (0,oo); 

(b) The integral identity 

r r 

u(Atp+\x\- 2 (p t )dxdt = (2.2) 



holds true for any test function ip £ C^°(R N x (0, oo)), and u(0) = uq. 
We have the following: 

Theorem 2.2. Let u$ £ L\{R N ), uq nonnegative. Then there exists a unique weak solution u to the 
Cauchy problem (|2.1j) such that 

\x\~ 2 u t , Au £ L loc (<2*), \x\~ 2 u t = An a. e. in Q*, 

w/iere = R w x (0, oo) \ {(0, t) : t > 0}. 
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This type of solution is usually referred as a strong solution (see [20]). The proof will adapt ideas 
from the proof of existence in [TOj Section 6] , where well-posedness is proved for (jl.4p with 7 > 2 and 
m > 1, thus at some points we will be rather sketchy. 

Proof of Theorem \2.2l Uniqueness. This is based on the following result which is also interesting by 
itself. 

Proposition 2.3 (L^-Contraction principle). Let u\, 112 be two strong solutions of Eq. (|l.ip . For 

< t\ < t2 we have 

WW.,*) -«,(«,<,)],*< / m-mm.) (2.3) 



where [g]+ represents the positive part. 

Proof. We follow the same ideas as in [201 Proposition 9.1]. Let p £ C^QR) n L°°(R) be such that 
p(s) = for s < 0, p'(s) > for s > and < p < 1. Let j be the primitive of p such that j(0) = 0. 
The idea is to choose p as an approximation of the function 



signer ) 



1, ifr>0, 
0, if r < 0, 



thus j will approximate the positive part function. Consider also a cutoff function £ n constructed in 
the following way: let £ € C$°(R N ), < £ < 1, £o0) = for |z| > 2, f (x) = 1 for |x| < 1, and 
define £, n (x) '■= Co(x/n). Subtracting the equations satisfied by U2 and u±, then multiplying by the test 
function p{u\ — U2)£, n , we obtain 



rt2 r m r 

/ / |x| _2 (ui - U2)tp{u\ - U2)Cndxdt = / / A(ui - U2)p(ui - u 2 )£„ dx dt. (2.4) 
Jh Jr n Jt x Jr n 

By approximation of u\ — U2 as in the proof of j20i Proposition 9.1] and integration by parts in the 
right-hand side of (|2.4p . we obtain, with ^ = £ n : 



[•12 r rt2 r 

I j \ x \~ 2 { u i — U2)tp{u\ — U2)i dx dt < — I I p{u\ — «2)V(tii — U2) • V£ dx dt 

Jt! Jr n Jti iR" 

Vj{u\ — U2) ■ V£ dx dt 

rt2 r rt2 r 

= / j (u\ — U2) A£ dx dt < / / \u\ — U2I I A£| dx dt. 

Jt! it* Jt! Jr N 

We let now p — > sign^ and integrate to get 

/ \x\~ 2 [ui(x,t 2 ) - U2(x,t 2 )} + £, n dx < / \x\' 2 [m(x,ti) - u 2 (x,ti)] + £ n dx 

JR N JR N 

+ ||A^ n ||oo / / \U1-U2\dxdt. 



Jt! J\x\>n 

Letting now n — > 00 and taking into account that we deal with solutions in L , we obtain the result. 
□ 

The uniqueness follows obviously from the contraction principle. Moreover, we have the following: 

Corollary 2.4 (Comparison principle). If u\, U2 are solutions of Eq. (jl.ip such that their initial data 
satisfy uq^i < uq^, then u\ < U2 in R N x (0,oo). 
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Proof of Corollary \2.4\ Suppose that Uo,i, U02 £ L>2(M. N ) satisfy Uo,i < ^0,2 and in, 112 are the corre- 
sponding solutions. Then, for any t > 0, we have 

/ \x\~ 2 [ui(x,t) — U2(x,t)] + dx < / |x| -2 [ito 1(3;) — ^0 2(2;)] 1 dx = 0, 
Jr n Jr n 

hence u\(x,t) < U2(x,t). □ 




Existence. This part of the proof is more involved and will be divided, for a better comprehension, 
into several steps. 

Step 1. Compactly supported data. In a first step, we prove existence of a strong solution for 
compactly supported initial data outside the origin. For 0<r<i?<oo, we denote an annulus by 
B r . jR := B R \ B r . Given uq £ Cq°(M n \ {0}), consider the approximating Dirichlet problem 

= An, in Q T)R := B r>R x (0, 00), 

■■ Uo(x), in B r>R , (2.5) 
0, in dB T)R x (0,oo). 

The problem (j2.5|) . being nor singular at x = neither degenerate at infinity, admits a unique solution 
u T)R as shown in [H [5]. Moreover, by simple comparison in the smaller annulus, it is easily seen that 
u r2 ,R 2 > u ri,Ri if r 2 < r i and i?2 > It remains to show that they are uniformly bounded. To this 
end, recall that uq 6 Cq°(R \ {0}) and consider r > small enough and i? > large enough such that 
suppiio C B r ^ R . Recalling the definition of F(x,t) given in (|1.8p . for some r > sufficiently small and 
K > large, we have 

ito(^) < KF(x,t) in i? ri R. 

Since on the lateral boundary u r R (x,t) = < KF(x,t + r) for any x such that either |x| = r or 
|x| = i?, by standard comparison we obtain the following important universal bound 

u rtR (x,t) < KF(x,t + r), in B r>R x (0, 00). (2.6) 

This is an uniform bound which does not depend on r, R. Thus, the following limit 

u(x,t) = \im lim u rR (x,t), (2-7) 

r— >0 _R->oo ' 

is well-defined and gives rise to a weak solution to (jl.lj) . as it is easy to check with the definition. In 
this way, we have solved the problem for no £ Cq°(M. n \ {0}). Moreover, let us notice here that any 
solution obtained through this approximation process also satisfies the universal bound (|2.6|) in Q*. 

Step 2. Preliminary estimates. We still want to prove that our weak solutions are in fact strong, 
then pass to the general case by some approximation process. To this end, we establish two estimates 
that will imply further regularity on the solutions. They are gathered in the following 

Lemma 2.5. Let u be a weak solution to dTTTJ) such that u G Cfi°(R N \ {0}). Let K CR N \ {0} be a 

compact set and < r < T < 00. Then there exists a positive constant C > depending only on N , 
dist({0}, if) and r such that 

T 

Vu\ 2 dxdt<C (2. 

K 

and 

pT p 

- 2 \u t \ 2 dxdt <C. (2.9) 

K 
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Proof. The proof adapts the ones in [10\ Theorem 3.6, Theorem 3.8], but we give a complete version 
of it for the sake of completeness and simplicity. Let K be a compact neighborhood of K such that 
does not belong to K and consider a cutoff function 77 £ C^°(M. N ) such that 

j] = 1 in K, i] = in \ K, < r) < 1 in R N . 

We multiply (jl.ip by urj 2 and integrate by parts to get 



/ \x\ 2 — 7] 2 dx = / uAurj 2 dx = — I V 'u ■ V '{urj 2 ) dx 
J~K 2 Jt^ Jt^ 



|Vti| 2 77 2 dx — 2 I ur/SJu ■ Vr/ dx. 



d f , ,_ 9 ti 2 

Integrating on [r, T] and using Young's inequality, we further obtain 

- / \x\~ 2 u 2 (x, T)r] 2 dx + / / | Vu\ 2 r/ 2 dx dt = - / \x\~ 2 u 2 (x, t)t] 2 dx — 2 / / ur/V u ■ V r] dx dt 
2 Jk~ Jt Jk: 2 J-^ J T 



< — [_\x\ 2 u 2 (x,t)t] 2 dx + 
2 J~K 



- f [_\Vu\ 2 r] 2 dxdt + 2 [ [_u 2 \Vr]\ 2 dx dt, 
2 J T Jk Jt Jk 



whence estimate (|2.8|) follows from the uniform bound (|2.6|) . which implies readily an uniform bound 
for the terms in the right-hand side. Consequently, 



I [ \Vu\ 2 dxdt< [ [_\Vu\ 2 r! 2 dxdt<C 
Jt Jk Jt Jk 



as desired. 

In order to prove the second estimate, we multiply (jl.ip this time by utrj, then integrate in space 
and time. We obtain 



whence 



I |-2| |2 i 

\x\ \ut\ rjax = - 

K JK 



rT 

-2|„. |2 



Vu ■ V(utrj) dx = —-— / \Vu\ 2 rjdx, 
2 dt Jj( 



x\ \ut\ dxdt< / / \x\ \ut\ rjdxdt 

Jt JW 
1 



,Vii(:E,T)| rjdx 

2 JW 

the last step coming from the first estimate (|2.8p . □ 



- [_\Vu(x,T)\ 2 i]dx < C, 
2 Jk 



In particular, we deduce from Lemma 12.51 that \x\ l ut = \x \Au £ L 2 oc (R N x (0,oo)), which at its 
turn shows that \x\~ 2 u t , Au G Lj oc (Q*). Hence u is a strong solution. 

Step 3. General data. For general data no £ L\(U. N ), we use the contraction principle (|2.3p . Let iio,n 
be a sequence of functions in C^°(R N \ {0}) such that UQ^ n — > uq in L\{R. N ) and u n the corresponding 
solution with initial datum uo, n - By (|2.3|) . the sequence u n is uniformly Cauchy in C([0, oo); L\{W N )), 
hence it converges to a limit u £ C([0, oo); L^IR^)) which is a weak solution (passing to the limit in the 
identity (|2.2p is immediate). Moreover, we deduce from (|2.6p that u is uniformly bounded on compact 
subsets of Q*. On the other hand, we apply Lemma [2.51 for u n and we obtain the estimates (|2.8p . (|2.9|) 
locally uniformly with respect to n in Q*. It follows first than the limit u satisfies the bound (|2.6p . 
then (|2.8p , and finally, by redoing for u the last step in the proof of Lemma I2.5|. we deduce that u also 
satisfies (|2.9p . Hence u is a strong solution. □ 
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3 Radially symmetric solutions. The transformations 

We restrict ourselves first to radially symmetric solutions u(r,t) = u(\x\,t), r = \x\, of Eq. (jl.ip . In 
this case, we introduce the following change of variable 

u{\x\,t) = v(y,t), y = log\x\ + (N-2)t. (3.1) 

By a simple calculation, we notice that 

u t (r,t) = (N - 2)v y (y,t) +v t (y,t), u r (r,t) = -v y (y,t), u rr (r,t) = \(v yy (y,t) -v y (y,t)), 

hence, by replacing in (jl.ip . we obtain that v is a solution to the Cauchy problem for the one-dimensional 
heat equation 

v t (y,t) = v yy (y,t), v (y) = u {r). (3.2) 

for (y, t) 6 R x (0, oo). We observe that x = is mapped through (|3.ip into y — > — oo. We thus can use 
the previous transformation (|3.1|) in order to prove Theorems II. II and II .21 But before doing that, let us 
also notice that, in dimensions N > 3 and also N = 1, there exists another transformation, appearing 
in a more general version in [161 Subsection 1.3.4], 

N — 2 

u(r,t) = e z ~ l w(z,t), z = — log|x|, (3.3) 

again mapping radially symmetric solutions u(r,t) of Eq. (jl.ip into solutions w(z,t) of (|3.2p . but with 
a different connection between initial data: 

«o(|x|) = e z w (z) = \x\^ N ^ 2 w (z). 

The change of variable and function (j3.3|) acts as an inversion for N > 3: it maps x = to y — > +oo 
and \x\ —> +oo to y — > — oo, and it acts as a direct mapping only for N = 1. In the sequel, we will use 
mainly transformation (|3.ip due to its simplicity, but at some points (|3.3p will appear too. 



4 Radially symmetric solutions. Asymptotic convergence 

In this section we prove Theorems 11.11 and 11.21 for radially symmetric solutions and we deduce some 
more interesting remarks about the quite unexpected asymptotic behavior for Eq. (jl.ip . 

Proof of Theorem li.il Let u be a solution for (jl.ip with initial datum satisfying (jl.5p and (jl.6p . By 
transformation (|3.ip . we arrive to a solution v of (|3.2|) with initial datum vq such that 



and 



lim v (y) = lim v (y) = 

y— >-oo y—t—oo 



M V0 := I v (y)dy= j dr 



-oo 

oo 



r 



(4.1) 



"^oWr^ 1 dr = — [ |a;r%,(|x|) dx = ^ < oo, 

where, as usual, oo\ is the area of the unit sphere of R . In fact, we notice that the same equality (j4.ip 
holds true taken at any time t > instead of t = 0. From the standard mass conservation property for 
the heat equation, we deduce that the quantity f RN \x\~ N u(x, t) dx is conserved. 
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Due to well-known results in the theory of the heat equation, we find 

1 



hm *V2 



t— >oo 



v(y,t)-M, 



V<3 



0. 



uniformly for We translate (|4.2p in terms of u to get 



lim t l/2 



t—>oo 



u(\x\,t) 



M,, 



1 



00l V47Tt 



cxp 



(log|x| + (N - 2)tf 
It 



(4.2) 



(4.3) 



that is (fTT7|) . □ 

We prove now the asymptotic convergence for the case when «o(0) = K > 0. 

Proof of Theorem li.ffi (a) Let u be a solution for (|l.ip with initial datum satisfying uo(0) = K > 0, 
(|1.5p and (|1 . 13[) . By the transformation ()3.ip . we arrive to a solution v of (|3.2p whose initial datum 
satisfies 

K > 0. 



lim u (y) = 0, 

y— >oo 



lim ■uo(y) 

j/-s— oo 



(4.4) 



Moreover, we notice that 

\KH(y) - v (y)\ dy 



\K - v Q (y)\ dy + 



f 

Jo 



1 \K-u (r) 



dr + 



My) I d y 



Un(r) , 1 

dr < — Jj < oo, 



(4.5) 



ff(y) 



due to the conditions on no in (|1.13p . where we have denoted by H the complementary Heaviside 
function 

1, fory<0, 
0, for y > 0. 

We now use the following 

Lemma 4.1. Let v\, t>2 be two solutions of the one- dimensional heat equation with initial data satisfying 

roc 

\v 2 (y,0) - Vl (y, 0)\ dy <oo. (4.6) 



Then we have 

\\v 1 (t)-v 2 (t)\\ 00 =0{t~ 1 / 2 ), ast^oo. (4.7) 

The proof of Lemma [4.11 is immediate, since w(y,t) := V2(y,t) — v±(y,t) is a solution for (|3.2p with 
integrable initial datum, thus has the desired order as t — > oo by standard results for the heat equation. 

We apply now Lemma 14 . 1 1 for the solutions v(y,t) and Ker£c(y,t)/2 of (|3.2p . Their initial data are 
vo(y) and KH{y) in the previous notations, thus they satisfy the condition (|4.6p . due to fj4.5|) and the 
fact that I\ < oo. Thus, we get that 



v (t) - ^erfc(t) 



0(t 



-l/2s 



as t 



oo, 



whence, undoing the transformation (|3.ip and coming back to the initial variables, we deduce (jl.lip . 

(b) Let u be a solution for (jl.ip with initial datum satisfying no(0) = K > 0, (|1.5p . (|1.13j) and 
(|1.14p . By the transformation ()3.ip . we arrive to a solution u of (j3.2|) with initial datum wo satisfying 
(|4.4p . We make a further change by letting w(y,t) := K — v(y,t), hence w is a solution to (j3.2j) with 

lim wo(y) = K > 0, lim w${y) = 0. 

y— >oo y— oo 

Let ip(y,t) := w y (y,t) = —v y (y,t). Then again V is a solution for (|3.2p . and we want to apply for ^ 
the following result: 
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Lemma 4.2. Let ip be a solution for Eq. (j3.2|) . whose initial datum ipo(y) := V'Cj/jO) satisfies the 
following conditions: 



:= / My) dy£ (0,oo), ^(V) : 

J— oo 

Then, the following two results hold true: 

U(y,t)-MWG(y,t)\\ P = 0(t- 2+1 / 2 n, 
where G(y,t) is the standard Gaussian 

G(y,t) 



'A-nt 



\y 3 ^o(y)\ d y < °°- 



(4>(s,t) - M(i>)G(s,t))ds 



Oit 



-3/2^ 



(4. 



(4.9) 



Lemma 14.21 is proved in |14l Theorem 2]. There, it is asked that tpo > 0, but it is easy to check that 
the result holds for any ip bounded. We have to check that tp satisfies the condition (|4.8p . We have: 



M(y>) 



V'o(y) dy = lim (w Q (y) - w Q (-y)) = K G (0, oo), 

i » 



and 



\y\{yM 



— oo 
oo 



I (log r) 3 ru r (r, 0)1 — 



/ |(logr) 3 n ,T-(^)| dr = — I 
JO w l JM 



l(10g ]l )3| |V^(|x|)|(ix<oo, 



due to the condition (|1.14p . Thus, by Lemma 14.21 we obtain the asymptotic convergence for ip as in 
3]), with M(ip) = K. We keep the second result in (|4.9p and transform it to get 



0(t~ 3 / 2 ), 



w(y,t) -K I G(s,t)ds 

J — oo 

or equivalently, recalling that w(y,t) = K — v(y,t), 

POO 

v(y,t) - K j G{s,t)ds 
Jy 

Noticing that 

"OO 

G(s,t)ds 



0(t~ 3 / 2 ). 



(4.10) 



1 



'Aitt 



e~ s2 / 4t ds 



2Vi 



2y/7Tt Jy/(2Vt) 



~ z2 dz 



-erfc 



2y/i 



and finally undoing transformation (|3.1|) to get back to the initial variables, we arrive to (|1.15p as 
desired. □ 

Remarks, (i) Decay rates as t — > oo. For x 7^ fixed, the behavior as t — > 00 of E(x,t), F(x,t) is 
different with respect to the dimension N. Indeed, when N > 3, both profiles decay to as t — > 00 
with a rate 



1 



exp 



(TV -2) 



as it is easy to check from the explicit formulas. Since the two special solutions make sense also for 
N = 1 and N = 2, we analyze the decay rate also in these cases. For N = 2, we have lim E(x,t) = 1 

t— >oo 
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and F(x,t) = 0(1/ y/t) as t — > oo, and for N = 1 we have lim E(x,t) = 2 and F(x,t) behaves in the 

t— >oo 

same way as for N > 3. 

(ii) Self-map. Consider Eq. (jl.ip in dimensions N, N respectively. Since both are mapped to solutions 
of (|3.2p through their correspondent transformation (|3.ip . we have 

y = log \x\ + (N- 2)t = log |x| + (N- 2)t, 

whence we obtain the mapping between radially symmetric solutions in the two dimensions 

u(x,t) = u(x,t), \x\ = |x|e (F_Ar)t . (4.11) 

The transformation (|4.1ip is a self-map of Eq. (jl.ip between dimensions N and iV. 

(hi) Hotspots. Inner and outer regimes. Since F(0,t) = and F(x,t) — > as t — > oo, it is natural 
to think about the evolution of the maximum points of F at time t (called hotspots). At this point we 
use the second transformation ()3.3p . in order to get, in the new variables (z,t), 

1 (JV-2) 2 , z 2 

F(x,t) = F(z,t) = -=e ~ l e WW*. 

V^TTt 

We omit the intermediate calculations which are straightforward. Notice that the maximum lies at 
zq = (N — 2) 2 t/2, and at that point, F(x,t) = l/V 47rf. Hence, the hotspots evolve with a decay rate 
which is smaller than the "linear" decay rate as t — > oo. This is a typical feature of the presence of two 
different regimes of convergence, one for \x\ small (the so-called inner regime) and another one for \x\ 
large (the so-called outer regime). This phenomenon exists also in literature for the porous medium 
equation in domain with holes in dimension N > 3 [T] or for the p-Laplacian equation in dimension 
N > p [9]. The difference between the behavior close to x = and close to infinity will become more 
clear in the case m > 1 that will be studied in the forthcoming paper [8]. 

5 Asymptotic convergence for general solutions 

We are now ready to prove Theorem 11.21 in its maximal generality. 

Proof of Theorem \1.2l Let u be a solution of Eq. (jl.ip with initial datum uq satisfying (|1.5p and (jl.lOp . 
with uo(0) = K > 0. We construct the sub- and supersolutions u_(x,t), u + (x,t) that are the (radially 
symmetric) solutions of Eq. (jl.ip having initial data 

u_.o( r ) := min{iio(x) : \x\ = r}, u+o(r) := max{«o(x) : \x\ = r}, (5-1) 

which are obviously continuous and bounded, and both satisfy (jl.lOp . thus also (jl,13p . Thus, by the 
comparison principle, Corollary 12. 4j we have 

u-(x,t) < u(x,t) < u + (x,t), (x,t) G M N x (0,oo). 

From part (a) of Theorem 1 1.31 we know that (|1.11|) holds true for U- and u + . Since u_ ) o(0) = n +j o(0) = 
K, it follows that 



K 

u(x,t) - —E(x,t) 



< max 



K 

u-{x,t) - —E(x,t) 



K 

u+{x,t) - —E(x,t) 



oir 1 ' 2 ) 



whence u satisfies (jl.lip . □ 

At this point we also prove Proposition 11.41 
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Proof of Proposition [L4\ (a) Let u be a radially symmetric solution with no = M5q. Then, by (|3.1D . 
its transformed v will be a solution to (j3.2j) . whose initial data will be the zero function, since x = 
is mapped into y — > — oo. But in our class of solutions, this solution of (|3.2p is the zero function, by 
standard uniqueness results Theorem 7, p. 58], which is a contradiction. 

(b) Assume in a first step that u is radially symmetric. We map u(r,t) into v(y,t) by (|3.1|) . Then 
vq = in (— oo,logro). Due to standard heat equation theory [6l Section 2.3.3 (a)], v(y,t) > for any 
t > 0, y £ (— oo,logro), but it remains true that lim v(y,t) = 0, for any t > 0. Coming back to it, 

y— >— oo 

we reach the conclusion. 

Let now u be a non-radial strong solution. We define the radially symmetric functions u_ , u + as in 
(|5.ip . Since tt_(r, £) < u(x,t), it follows that u(x,t) > for any x / 0, i > 0. Since u(x,t) < u + (r, i), 
it follows that it(0, t) = 0, for any i > 0. □ 



U(x,t) 



6 Extensions, comments and open problems 

1. Dimension iV = 1. In this case, the transformation (|3.ip applies to any solution. We thus notice 
that the effect of the singularity at x = is just disconnecting the real line. More precisely, let u be a 
continuous solution to (jl.lj) posed in dimension N = 1 and such that u(0, t) = for all i > 0. Define 

/ ,v f u(x,t), if x > 0, / .% f u(x,t), if x < 0, 

U+(x ' t):= \0, ifx<0, M:M):=| 0) ifx > Q) 

and notice that, applying the transformation (|3.ip to both it + and u_, we obtain two different solutions 

v + (y, t) := u+(x, t), y = log x — t, x > 0, and v-(z, t) = U-(x, t), z = log(— x) — t, x < 0, 

which will asymptotically converge towards two different Gaussian profiles with different masses M + 
and M_ as indicated in the proof of Theorem 11.11 Coming back to initial variables, the profile of u 
will be 

' M + F(x,t), if x > 0, 
M-F(x,t), ifx < 0. 

In particular, we give a simple explicit example of solution having two branches of the type indicated 
above: 

P/ w l._/ a %')i forx<0, . . 

l[,, -~\(l-a)F(x,t), forx>0, [ 1 

for any a £ (0, 1), a ^ 1/2. We leave the details of the proofs of the previous statements to the reader, 
along the lines of the proof of Theorem 11.11 For initial data with uq(0) = K > 0, Theorem 11.21 still 
applies in this case. In order to show the difference of behavior with respect to N > 3, we plot in 
Figure [2] the two tipical profiles for dimension N = 1. 

2. Improvements and open problems related to Theorems 11.21 and 11.31 Besides the conditions 
in Theorem 11.31 in the case uq(0) = K, one can improve the rate of convergence to the asymptotic 
profile E(x,t) by using, for example, results from the paper [21]. We leave these improvements to 
the reader, along the lines of the proof of Theorem 11.31 On the other hand, an open problem arises 
naturally from the results in Theorems 11.21 and 11.31 That is, can one prove (jl.lip for general (that 
means, not necessarily radially symmetric) solutions, only asking the condition (|1.13p to hold? 

3. Asymptotic profiles of non-radially symmetric solutions with uq(0) = 0. Theorem 1 1 . 1 1 and 
the counterexamples following it raise a very natural, but also in our opinion very difficult problem 
to classify the asymptotic profiles for general solutions of Eq. (jl.ip when uq(0) = 0. Since the 
counterexamples are constructed in a very specific way, a first partial question that could be raised is: 
are there all the profiles weighted combinations of the radial one? For dimension N = 1 we give an 
answer in the first comment of this section, but the problem remains open for the rest of dimensions. 
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Figure 2: Profiles E and F for dimension N = 1 at various times. 



4. General densities. As we have said in the Introduction, there exists an increasing interest for the 
more general problem (|1.3p . This case has not been yet studied in detail with g(x) ~ \x\~ 2 as |x| — > oo 
and m = 1, but it has been strongly studied for m > 1, [TSIEHHU] and references therein. The authors 
consider there densities g which are regular near x = and show that the asymptotic behavior is given 
by the fundamental solution to (jl .4|) . It might be interesting to raise also the problem of considering 
general densities of various other forms, for example 

g(x) ~ \x\~ 2 as \x\ — > 0, lim g(x) = C G (0, oo), 

that is, preserving the other property of the pure power density (the singularity at x = 0) and renounc- 
ing to the first one (the behavior as \x\ — > oo). 
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